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In the factorized picture of semi-inclusive deep inelastic scattering the naive time reversal- 
odd parton distributions exist by virtue of the gauge link which is intrinsic to their 
definition. The link structure describes initial/final-state interactions of the active par- 
ton due to soft gluon exchanges with the target remnant. Though these interactions are 
non-perturbative, calculations of final-state interaction have been performed in a per- 
turbative one-gluon approximation. We include higher-order contributions by applying 
non-perturbative eikonal methods to calculate the Boer-Mulders function of the pion. 
Using this framework we explore under what conditions the Boer Mulders function can 
be described in terms of factorization of final state interactions and a spatial distortion. 

Keywords: transverse momentum dependent parton distribution function (TMD); semi- 
inclusive deeply inelastic scattering (SIDIS); single spin asymmetry (SSA) 
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1. T-odd PDFs, Gluonic Poles and The Lensing Function 

Over the past two decades the transverse partonic structure of hadrons has been 
the subject of a great deal of theoretical and experimental study. Central to these 
investigations are the early observations of large transverse single spin asymmetries 
(TSSAs) in inclusive hadron production from proton-proton scattering over a wide 
range of beam energies. Recently TSSAs have been observed in lepton-hadron 
semi-inclusive deep inelastic scattering (SIDIS)^"^ as well as in inclusive produc- 
tion of pseudo-scalar mesons from in proton-proton collisions. ^'^"^'^ While the naive 
parton model predicts that transverse polarization effects are trvial in the helicity 
limit/^ Efremov and Teryaev demonstrated^^' that soft gluonic and ferniionic 
poles contribute to multiparton correlation functions resulting in non-trivial twist- 
three transverse polarization effects in this limit. ^^'^^ In addition theoretical work 
on transversity^^"^^ indicates that transverse polarization effects can appear at lead- 
ing twist. Two explanations to account for TSSAs in QCD have emerged which are 
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based on the twist-three^^' and twist-two^^"^^ approaches. Recently, a coherent 
picture has emerged which describes TSSAs in a kinematic regime where the two 
approaches arc expected to have a common description. ^^"^^ 

In the factorized picture of semi-inclusive deep inelastic scattering^'*' ^° at small 
transverse momenta Pt ^ kr « \/Q^ the Sivers effect describes a transverse 
target spin-S'x asymmetry through the "naive" T-odd structure, Af{x,kT) ~ 
St ■ (P ><- kT)f^{x,k'^).^^'^^ For an mipolarizcd target with transversely polar- 
ized quarks-ST, the Boer-Mulders function-^^ is Ah{x, kx) ^ st ■ {P ^ kT)hi{x, k^)- 
Many studies have been performed to model the T-odd PDFs in terms of the FSIs 
where soft gluon rescattering effects are approximated by perturbative one-gluon ex- 
changed^"''* We improve this approximation by applying non-perturbative eikonal 
methods to calculate higher-order gluonic contributions from the gauge link in the 
spectator framework. In the context of these higher order contributions we per- 
form a quantitative study of approximate relations between TMDs and GPDs. In 
particular, we explore under what conditions the T-odd PDFs can be described via 
factorization of FSI and spatial distortion of impact parameter space PDFs.*^ While 
such relations are fulfilled from lowest order contributions in field-theoretical spec- 
tator models'**' a model-independent analysis of generalized parton correlation 
functions (GPCFs)'*^ indicates that the Sivers function and the helicity flip GPD E 
are projected from independent GPCFs. A similar result holds for the Boer-Mulders 
function for a spin zero target.'*^ From phenomenology, however it essentially un- 
known whether the proposed factorization is a good approximation. Here we focus 
on the transverse structure of the pion through the impact parameter GPD and 
the Boer Mulders function for which very little know. Recent lattice calculations 
indicate that the spatial asymmetry of transversely polarized quarks in the pion is 
quite similar in magnitude to that of quarks in the nucleon.** 

The field-theoretical definition of transverse-momentum dependent (TMD) par- 
ton distributions in terms of hadronic matrix elements of quark operators for spin- 
1/2 hadron with momentum P and spin S was presented in Refs.^*'^^''*^ It is 
straightforward to obtain the TMDs for a spin-0 hadron from that. One encoun- 
ters two leading twist TMDs for a pion, the distribution for unpolarized quarks /i , 
and the distribution of transversely polarized quarks hj^, the Boer-Mulders func- 
tion. Adopting the infinite-momentum frame where the hadron moves relativistically 
along the positive 2:-axis such that the target momentum P has a large plus compo- 
nent P+ and no transverse component the Boer-Mulders function, defined in SIDIS 
is 

The light-like vector n represents a specific direction on the light-cone n'^ = (1, 0, 0) 
where we define the light cone components of a 4- vector = l/\/2{aP ± a^), 
= (a~, a^, a"*"), [x ; y] denotes a gauge link operator connecting the two locations 
X and y. We work in a covariant gauge where a transverse gauge link at light-cone 



December 30, 2009 1:28 



Gamberg 



Final state interactions and the transverse structure of the pion 3 

infinity is negligible. The gauge link in ^ is interpreted physically as FSIs of the 
active quark with the target remnants^^'^*^ and is necessary for "naive" time-reversal 
odd (T-odd) TMDs^^^ 25,31 exist.^o The Boer-Mulders function appears in the 
factorized description of semi- inclusive processes such as SIDIS2**"26,30, 51-56 terms 
of the first fcT-moment, 2m'^hf^'^\x) = J (P kx hj^ {x , k^) . Transforming the two 
pion states in Eq. ([1]) into a mixed coordinate-momentum representation specified 
by the impact parameter 5t*^'^^ results in an impact parameter representation for 
the gluonic pole matrix element^^ {kT){x) = rriT^hf^^^x), 

(krHx) = J cmT^^e"^^'''{P+,0T\q{zi)[zi;z2]r{z2y+q{z2)\P+,0T). (1) 

Here, the impact parameter is hidden in the arguments of the quark fields, 
zf^^2 — T^n-^ + and 6^ = {0,b}p,b^,0). The operator /* originates from the 
time-reversal behavior of the FSIs written in terms of the gauge link in H]) and the 
field strength tensor F^" , 

2/^(^2) - I dy- [Z2 ; y] gF+\y) [y ; z^], (2) 

with y^^ = y-n^" -I- 6^. 

Turning our attention to GPDs of a pion, they are represented by an off-diagonal 
matrix element of a quark-quark operator defined on the light-cone. One en- 
counters two leading twist GPDs for a pion, a chirally-evcn GPD Fi and the chiral 
odd GPDs Hi ^"^ We use the symmetric conventions for the kinematics for GPDs,^^ 
P = i(p + and A = p' - p where A+ = -2^P+, and t = A^. The impact 
parameter GPDs are obtained from the ordinary GPDs via a Fourier-transform of 
the transverse momentum transfer A^ at zero skewness ^ = 0. The chirally-odd 
impact parameter GPD Wl is expressed as 

/rl7~ 4- - ^ W r) ^ 

2(27r) db^ 

Ti.1 describes how transversely polarized quarks are distributed in a plane transverse 
to the direction of motion. This distribution represents transverse space distortion 
due to spin-orbit correlations.^*' ^^'^^ A comparison of the first moment of the Boer 
Mulders function ([T]) and the impact parameter GPD TCi reveals that they differ by 
the operator /* which represents the FSIs. In various model calculations, the FSIs 
are approximated such that^^'**^' 57,63 |.]^g effects of a distortion of the transverse 
space parton distribution and the FSIs factorize resulting in the quantitative relation 

m^elfht'-'\x) [ d'bTlHx,bTf-^ -i-ni{x,bl), (4) 
J iTiiT db^ 

where I is the so-called "quantum chromodynamic lensing function" .^^ This factor- 
ization doesn't hold in general. ^^'^"^ For example, this relation breaks down when 
the quark fields and the operator / in ([1]) "interact" via quantum fluctuations (be- 
cause they are interacting Heisenberg operators). However, it unknown if ([4| is a 
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Fig. 1. The amplitude W including FSIs between re-seattered eikonalized quark and antiquark. The 
FSIs are deseribed by a non-perturbative scattering amplitude M that is calculated in a generalized 
ladder approximation. Gluon interactions as shown in the second diagram are not taken into account 
(see text). 

good phenomenological approximation. We estimate the size of the lensing function 
using non-perturbative eikonal methods^^'^^ to calculate higher-order gluon contri- 
butions from the gauge link. Up till now the relation ([4]) was used to predict the sign 
of T-odd TMDs in conjunction with numbers for the u- and d-quark contributions 
to the anomalous magnetic moment of the nucleon and the assumption that final 
state interactions are attractive. 



2. TMD - GPD Relation for a Pion the Spectator Framework 

We focus our attention in the following on a pion in a valence quark-type config- 
uration that one would expect for relatively large Bjorken x. Thus, we only take 
valence quark wave functions into account and neglect higher Fock states. Assuming 
an antiquark spectator one can express the pion Boer-Mulders function ([1]) as 

where W°''^{P, k; a) = {P — k, a, 6\ [oon ; 0]"'' (0) \P) and a and 5 represent the 
helicity and color of the intermediate spectator antiquark. We model the matrix 
element for by a diagram shown in Fig. [1] where the final state interactions - 
generated by the gauge link in - are described by a non-perturbative amputated 
scattering amplitude (M)"^ with (3, a (7, S) color indices of incoming and outgoing 
quark (anti-quark), respectively. We neglect classes of gluon exchanges in the second 
diagram in Fig. [1] represented by the gluon rungs since they would be attributed to 
the "interaction" between the quark fields and the operator / in ([1]). They lead to 
terms which break the relation (j4]). We also neglect real gluon emission and (self)- 
interactions of quark and antiquark lines the second diagram in Fig. [T] since they 
represent radiative corrections of the GPD and are effectively modeled in terms 
of spectator masses and phenomenological vertex functions. The pion-quark vertex 
is governed by the interaction Lagrangian C — ~gTr/\^T%6"^q"j5T- and we 
allow the coupling constant to depend on the momentum of the active quark in 
order to account for the compositeness of the hadron and to suppress large quark 
virtualities. Phenomenological vertex functions in connection with spectator models 
have been used frequently in the literature.^^''*^'^^ Applying the Feynman rules it is 
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then possible to find an expression for the matrix element W from the first diagram 
in Fig. [2 We find 



where Ps = P — k is the spectator momentum. The first term in ([6]) represents 
the contribution without final state interactions while the second term corresponds 
to the first diagram in Fig. [T] We then express the FSIs through the amputated 
quark - anti-quark scattering amplitude M. Here both incoming quark and anti- 
quark are subject to the eikonal approximation (see, e.g. Rcf. 68 and references 
therein). While the active quark undergoes a natural eikonalization for a mass- 
less fermion since it represents the gauge link contribution, the eikonalization for a 
massive spectator fermion is a simplification that can be justified by the physical 
picture of partons in an infinite momentum frame. The eikonalization of a massive 
fermion can be traced back to the Nordsieck-Bloch approximation^^ which describes 
a highly energetic helicity conserving fermion undergoing multiple scattering un- 
dergoing very small momentum transfer. In this approximation the Dirac vertex 
structure, u{pi)"f'^u{p2) ^ Pfi/n^ where (pi + p2)/2 = p. For a massive anti-fermion 
one identifies the velocity = —p^/m, and the numerator of a fermion propagator 
becomes j(— ^ -I- m) — > i{—n ■ p -I- to). 

We proceed by performing a contour-integration of the light-cone loop- 
momentum, in Eq.Q. In doing so we only consider poles which originate from 
the denominators in fH]). This assumes that the scattering amplitude M does not 
contain poles in . This assumption is not necessarily true even for a one gluon 
exchange approximation. As shown in Refs. 70, 40, q^ poles appearing in M are re- 
lated to light-cone divergences that may be regulated by choosing a slightly off-light 
like vector n. Performing the contour integration under the described assumptions 
fixes the momentum q~ of the anti-quark in the loop in ^ to q^ — [q^ + rrij.) / 2q'^ . 
The eikonal propagator can be split into a real and imaginary part via the usual 
principle value prescription l/(x -I- iO) = V(l/x) — iTrS{x). It has been argued^^ 
that only the imaginary part contributes to the relation ^ as it forces the an- 
tiquark momentum q to be on the mass shell. Thus, the imaginary part of the 
eikonal propagator corresponds to a cut of the first diagram in Fig. [1] before the 
final-state interactions take place, whereas the real part given by the principle value 
is attributed to terms that break the relation ([4]). 

After performing these integrations we use ([6]) to calculate the pion Boer- 
Mulders function via ([5]) . The pion-quark-antiquark vertex is 




9. {{P ~ q?) {r-i + "ig)75(-^ + m,) {M)% {q, Ps)v{Ps,a) 



[n ■ {Ps -q) + iO] [{P - g)2 -m^+ iO] [g2 - to^ -|- iO] 



-I ^ 



(6) 



5'7r(fc^) = g-. 




{k^-ml)f(k^) 



(7) 



fc2 - A2 + ?:o 
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where the function / is a homogeneous function of the quark virtuaUty. We choose it 
to be a Gaussian exp[— |fc^ |] .""^ Inserting ([6]) into ([5]) yields the following expression 
for the Boer-Mulders function, 



{xrus + (1 - x)mq) ((1 - x)A'^y 



d?qT (PpT 
(27r)2 (27r)2 



-TIM' 



ciki 



(8) 



where 

J^[M^'^] ^ [kT + qr) ((2^)25"/55(2) (pV + kT) + {kr +Pt) 

(9) 

with A2(x) = xm^ — x{l — x)M^ + (1 — x)A^. This result already anticipates an 
eikonal form for the scattering amplitude M{x, kr, qr) M'^^^{\qT + krl) which we 
exploit to simplify the expression and show a relation to the chirally-odd GPD H^. 
A calculation for the GPD for an antiquark spectator can be found in Rcf. 47, 
which we generalize with the vertex function g-^{k'^) 



HI {x, 0,-^1) 



;\xms + (1 - x)mq) 



2n-l 



2(27r)3A2 



Z32 +A2(.t) 

2a2(d2 4.a2(x)) 



d^p I dz 



1 — 4z(l — Z)^ 3-— C0S2 03 



(10) 



where D^j. = i(l — xYA"^. Weighting ([5]) with a transverse quark vector kip and 
integrating both sides over kT we readily obtain the relation 

2 



mlh^^'\x) 



2(27r)2 



qT-Iix,qT)H^X,0 



gT 
1~X 



(11) 



The function X* can be expressed in terms of the real and imaginary part of the 
scattering amplitude M, 



Tix.qr) 



1 



d'^PT 



{2pT-qTy {m''''%U\PT\) 



NcJ (2^)2 
((2^)2r'3^(2)(^^ _ ^ I'sjjj^;^. 



07 



(Ipt - grl) 



(12) 



In order to derive the relation ^ one transforms Eq. pip into the impact parameter 
space via a Fourier transform. The lensing function in the impact parameter space 
then reads, 



X\x,bT) =i{l~x) 



(13) 
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3. The Lensing and Boer Mulders Function in Relativistic Eikonal 
Approximation 

In order to calculate the 2^2 scattering amplitude M (needed for (fT^ 'l we use 
functional methods to incorporate the color degrees of freedom in the eikonal limit 
when soft gauge bosons couple to highly energetic particles on the light cone. Here 
we summarize implementation of the color structure in the calculation of M while 
the details of the functional approach can be found in a forthcoming publication. 
Work in this direction was carried out in Refs. 71, 66, 72. After some functional 
manipulations of the scattering amplitude, M can be expressed in terms of quark- 
and antiquark-propagators that are linked together by soft colored gluons. The 
amplitude M reduces to a simple, gauge-invariant expression^^ 

{M'^-'X'^ (x, Wt + krl) - ^^^^ / d^^Te--^-('--+^^-) (14) 



(2^)^?-l \ 



In Eq. (HD) the A^ — 1 dimensional integrals over the color parameters results from 
auxiliary fields q;°(s) and u°(s) that were introduced in the functional formalism of 
Ref. 66 in order to decouple the gluon fields from the color matrices. The eikonal 
phase x(|2t|) in Eq- p4)) represents the arbitrary amount of soft gluon exchanges 
that are summed up into an exponential form and is expressed in terms of the gluon 
propagator in a covariant gauge, 

/C30 j'OO 
da I dPnf'n''Vf,^{z + an- /3n). (15) 
-OO J —OQ 

T) denotes the gluon propagator, and g the strong coupling. In this form the four- 
vector V is related to the complementary light cone vector n,v = —{{l—x)P~^/ms)n, 
with n-n — 1 and = 0. We evaluate the color integrals by deriving a power series 
representation for the color function 

After manipulating the exponentials in (|16p , rewriting the resulting factors as deriva- 
tives with respect to u, and performing an integration by parts we obtain the fol- 
lowing power series representation for /, 

fM = Y.7^ E - E E(*"-*""*"^"'^'-*"^"'"')./3' (17) 

n=l ^ ai = l a„ = l P„ 

where P„ represents the sum over all permutations of the set {1, ...,n}. If we had 
a direct ladder where gluons were not allowed to cross we would have only factors 
(f°i...f°"i""...i"i)Q/3 = CpSap with Cf = 2'Af ' ^^'^ could work in an Abelian 
theory with an effective replacement a — s- Cpcts for the fine-structure constant. 
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4 4 

Fig. 2. The function C[-j] of Eq. (I22t as a function of the ciitonal amplitude . Left: SU{2), we compare 
the numerical result computed by means of Eq. il7l up to the order n — 8 with the analytical result in 
Eq. ]24t . The numerical agrees with the analytical result up to ^ ~ 2. Right: SU{3)^ we compare the 
numerical results for the orders n —7,8. The results arc accurate up to ~ 1.5. 



Since we allow generalized ladders with crossed gluons we have to sum over all 
permutations in ()17p . and the simple replacement is not possible. In a large Nc 
expansion the crossed gluons diagrams would be suppressed such that the direct 
ladder represents the leading order in l/Nc- In an Abelian theory, the generating 
matrices t reduce to identity and since we have n! permutations of the set {1, n}, 
we recover the well-known Abelian result, 

n=l 

For Nc — 2, t"- — a"" /2 and we can calculate the integral analytically by means 
of the relation (e™"^)^^ = 6ap cos (^^^ + shi (^^^ ■ We obtain, 

C^'^ (f) = {^^' f " I " ^) + '^"f" + f l) ■ 

We also calculate numerically the lowest coefficients in the power series ([T71) . and 
they exactly agree with the coefficients in an expansion in x of the analytical result 
([TO)). This serves as a check of both numerical and analytical approaches. For = 3, 
due to difficulty of integrating over the Haar measure we use the power series ([T7|) 
to and obtain the following approximative color function which is valid if a = x/4 
is small, 

= Sc,p{-C2a^ + - cea'' ~ cga^ + ...), 
S[/fJ^"](a) = Sc,p{cia ~ c^a^ + c^a^ - cW + ■.■), (20) 

with the muiierical values ci = 5.333, C2 = 6.222, C3 = 3.951, C4 = 1.934, C5 = 0.680, 
cg = 0.198, C7 = 0.047, C8 = 0.00967. Working in coordinate space where we can 
express the lensing function directly in terms of the real and imaginary part of the 
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color function / which is itself a function of the eikonal phase x Eq. 
in a lensing function of the form, 

(1 - x) x' 



r{x,bT) 



2N^ \br\ 4 



■C 



with 



C 



'X' 




A. 





(Tr3[/])' (I) + -Tr (|) (|) 

m]) (f ) m])' (f ) 



5- 



results 



(21) 



(22) 



where x' denotes the first derivative with respect to \zt\, and (3[/])' and (5i[/])' 
are the first derivatives of the real and imaginary parts of the color function /. 
Inserting ([18]) into ((2T|) yields the lensing function in an Abelian ?7(l)-theory 



)(x, br) = (1 - ^)-^x'(^) 1 + cosx(^) . 
Likewise by using pop the lensing function in an SU (2)-theory is given by 

( 



T: 



5(7(2; 



(-.M = ^^::4^x'f3(i + cosJ)+^^ 



16167 



■ X fx ■ X 

sm — sm — 

4 V4 4 



(23) 



, (24) 



where X — X (^i^j ■ S'[/(3)-theory, we use the approximate color function / 

to calculate the lensing function as a function of the eikonal amplitude. In Fig. [2] the 
function C[j] is plotted versus j for various approximations. While the convergence 
of the power series seems to be better for SU{2) than in the SU{3) case where the 
numerical result calculated with eight coefficients agrees with the analytical result 
up to ^ ~ 2, we can trust the numerical result computed with eight coefficients up 
to f - 1.5 for S'[/(3). 

In order to numerically estimate the lensing function and in turn the Boer Mul- 
ders function we seek to utilize the infrared behavior of the gluon and the run- 
ning coupling in the non-perturbative regime where we infer that the soft gluon 
transverse momentum defines the scale at which the coupling is evaluated. These 
two quantities have been extensively studied in the infrared limit in the Dyson- 
Schwinger framework^^ and in lattice QCD.^** We use calculations of these quanti- 
ties from Dyson-Schwinger equations^^ where both as and 'D~^ are defined in the 
infrared limit (details can be found in a forthcoming publication) . This determines 
the eikonal phase and thus the lensing functions ((2T|) for a U{1), SU{2) and SU{3) 
color function. We plot the results in Fig. [3] for a color function for U{1), SU{2), 
SU{3). While we observe that all lensing functions are attractive and fall off at large 
transverse distances, they are very different in size at small distances. 

Using the eikonal model for the lensing function together with the spectator 
model for the GPD we present predictions of the relation ([4]) for the first moment 
of the pion Boer-Mulders function h^^^\ We fix the six free model parameters m^, 
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Fig. 3. Left: The lensing function J'(x, br) from Eq. ^ for (7(1), SU{2) and SU{3) for x = 0.2 
at a scale Agfjjj = 0.2 GeV. For comparison we also plot the perturbative result with an arbitrary 
value for the coupling, a = 0.3. Right: First moments of the pion Boer-Mulders function calculated 
by means of the relation to the chirally-odd GPD TiJ for a SU{3), 5(7(2), and U{1) gauge theory. 



rrig. A, A, g^^ and n in (jlOp that we need to determine by fitting to pion data . In 
order to do so we determine the chiral-even GPD in the spectator model by 
investigating different hmits of . When integrated over x, the GPD reduces to 
the pion form factor F"^* (Q^) — —F'^ (Q^)- An experimental fit of the Pion form 
factor to data is presented in Refs. 75, 76, and up to = 2.45 GeV^ a reasonable 
fit to the data is displayed by the monopole formula Ffit(Q^) = (1 + 1.85 
This procedure is expected to predict the t-dependcncc of the chirally-odd GPD 
reasonably well up to — 2.45 GeV^. In order to fix the a;-dependence of we 
fit the collinear limit F^ (x,0,0) to the valence quark distribution in a pion, v{x). 
A parameterization for this object was given for example by GRV in Ref. 77 at 
a scale /it^ = 2GeV^. Reasonable agreements of the Form Factor- F^{x,Q,0) with 
the data fits were obtained for the parameters rriq ~ 0.834 GeV, rus — 0.632 GeV, 
A = 0.067 GeV, A = 0.448 GeV, n = 0.971, = 3.604. With the predicted GPD 
Hi and the lensing function I*(a;, 6t) as input we can use the relation ([4]) to give 
a prediction for the valence contribution to the first fcx-moment of the pion Boer- 
Mulders function, which we can write as, 

roc Ci 

mlht^'\x) = 27: dbTbll{x,bT)Q^ni{x,hl)- (25) 

We present numerical results for xnn^ hi^^\x) shown in Fig.[3]for a ;7(1), SU{2) and 
SU{3) gauge theory.''^ It was argued in Ref. 57 that a negative sign of the lensing 
functions indicates attractive FSIs. We find that the lensing function is negative for 
the both the Abelian and non-Abelian gauge theories. The magnitude of the SU{3) 
result is about 0.01, while the SU{2) result and U (1) result are smaller. One observes 
a growth of the pion Boer-Mulders function with Nc which was also predicted by 
a model-independent large Nc analysis.^* So far the pion Boer-Mulders function is 
unknown but maybe determined from a future pion-proton Drell-Yan experiment 
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to be performed at COMPASS. Once a pion Boer-Mulders is extracted our analysis 
can be used to verify quantitatively GPD - TMD relations. An extraction of the 
other T-odd parton distribution, the proton Sivers function ,f^x^\ from SIDIS data 
measured at HERMES and COMPASS reveals an effect of the magnitude of about 
0.04, four times larger than our prediction. A similar calculation for the proton 
Sivers function will be reported elsewhere. 
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